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The pairing in neutron matter has been investigated in presence of realistic two- and three- 
nucleon interactions. We have adopted the Argonne tig' NN and the Urbana IX 3N potentials. 
Quantum Monte Carlo theory, specifically the Auxiliary Field Diffusion Monte Carlo method, and 
Correlated Basis Function theory are employed in order to get quantitative and reliable estimates of 
the gap. They both fully take into account the medium modifications due to the interaction induced 
correlations. The two methods are in good agreement up to the maximum gap density and both 
point to a slight reduction with respect to the standard BCS value. In fact, the maximum gap is 
about 2.5 MeV at k F ~ 0.8 fm" 1 in BCS and 2.3-2.4 MeV at k F ~ 0.6 fm" 1 in correlated matter. 
At higher densities the Quantum Monte Carlo gap becomes close to BCS. In general, the computed 
medium polarization effects are much smaller than those previously estimated within all theories. 
Truncations of Argonne t> 8 / to simpler forms give the same gaps in BCS, provided the truncated 
potentials have been refitted to the same NN data set. Differences among the models appear in the 
correlated theories, most of the reduction being attributable to the tensor force. The three-nucleon 
interaction provides an additional increase of the gap of about 0.35 MeV. 
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Superfluidity in highly asymmetrical nuclear matter, 
with a low concentration of protons, can lead to major 
consequences in modern astrophysics. 1 Sq pairing may 
occur in the low density neutron gas in the inner crust 
of neutron stars and in the proton phase in the interior 
of the star. Anisotropic 3 P2~ 3 F2 neutron pairing may 
also appear at higher interior densities. Understanding 
microscopically such phenomena as the cooling rate 0, 
and the post-glitch relaxation times 0, Q requires as 
accurate as possible knowledge of the properties of the 
phases (normal and superfluid) of nuclear matter. 

The role that tensor force and N-N correlations have 
on the superfluid properties of neutron matter of T = 
is not yet fully understood. In this paper we address 
this problem by performing ah initio calculations of the 
1 Sq BCS-type pairing in neutron matter, using (i) a re- 
alistic Hamiltonian and (ii) a trustworthy many-body 
theories. As far as the Hamiltonian is concerned, a num- 
ber of nucleon-nucleon interactions have been derived in 
the last decade, all of them fitting the large body (~ 
3000) of the available pp and np scattering data with a 



X 2 /datum ~ 1. The Argonne Via Nijemegen Q, and 
CD Bonn 7j interactions are members of this family of 
phase-equivalent potentials. The situation is much less 
established for the three-nucleon interaction (TNI) be- 
cause of the much smaller number of experimental data 
to be fitted as well as the larger uncertainty of the theo- 
retical framework on which its construction is based 

There has been a rapid advance in dealing with 
strongly interacting systems in their normal phase within 
modern many-body theories. Parallel to the more tra- 
ditional BHF and CBF theories there has been signifi- 
cant development of the Quantum Monte Carlo meth- 
ods (QMC) in dealing with spin-dependent Hamilto- 
nian. In the different versions of QMC the normal phase 
many-body Schrodinger equation is solved via stochas- 
tic sampling. For homogeneous systems, a finite number 
of particle is put inside a simulation box with periodic 
boundary conditions to simulate the infinite gas. Then 
either the spatial configurations only (Green's Function 
Monte Carlo, GFMC [l(J]) or both the spatial and spin- 
isospin ones, by introducing auxiliary fields (Auxiliary 
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Field Diffusion Monte Carlo, AFDMC PJ), are sam- 
pled. GFMC calculations are limited to a low number 
of particles in the box (N=14 in neutron matter) since 
it sums over all the 2 N spin states of the N interacting 
neutrons, whereas AFDMC does not have this limitation. 
As a consequence, GFMC has typicall y la rge finite size 
corrections. GFMC [HI and AFDMC [HI have been so 
far used to study pure neutron matter with comparable 
results. Actual calculations restrict the sampling within 
the nodal surface of some trial wave function, where the 
sampling guiding function has a constant sign (Path Con- 
straint approximation 01). An accurate choice of the 
nodal surface is essential in obtaining good quality re- 
sults 0. This is more difficult for AFDMC since the 
method samples the spin degrees of freedom, and has 
less inherent cancellation. 

Correlated Basis Functions (CBF) theory provides 
an alternative way of addressing interacting nuclear sys- 
tems. The strong nuclear interaction modifies the short 
range structure of the wave function (short range corre- 
lations, SRC) and introduces many-body contributions. 
Within CBF the SRC are introduced through a many- 
body correlation operator acting on a model function 
(Fermi gas for normal phase infinite systems or shell 
model wave function for finite nuclei). Since the NN 
potentials have important spin- and isospin-dependent 
components, the nuclear correlation operator is highly 
state dependent. An effective correlation operator has 
been shown to be: 



F 6 (1,2,..N) = S 



i<j=l,N 



where 



/e(y) = E /^(ryJOWfo') 

p=l,6 



(1) 



(2) 



with 0^ 2 ^(ij) = 1, aii) ■ a(j), S(ij) = 
(3r a (i)rp{j) - 5 a/} )a a (i)ap(j), and C>(*> -*>+ 3 ) (£j) = 
O^iij) x [r(i) • t(J)], S being a symmetrizer operator 
since the fe operators do not commute among each other. 
In pure T = 1 neutron matter, the isospin components 
may be actually embodied into the spin dependent ones. 
This choice of the correlation is consistent with the use 
of a momentum independent vq potential: 



= E^, 



)0<">(y) 



(3) 



p=l,6 



Modern potentials also contain momentum dependent 
(spin-orbit and so on) components. Accordingly, spin- 
orbit correlations have been used in CBF nuclear matter 
studies . The correlation is variationally fixed by min- 
imizing the ground state energy as computed via cluster 
expansion. 



The normal phase of neutron matter has been studied 
in details by both QMC and CBF. Very recently we have 
extended the two methods to deal with the 1 S'o superfluid 
pairing. A detailed account of the theories is in prepara- 
tion and will be given in subsequent publications |l7lll8|. 
We present here results obtained within the AFDMC and 
the CBF methods for the gap energy computed starting 
from realistic potentials (Argonne vi$ and its reductions 
plus the Urbana IX TNI) and compare with the pure 
BCS results. A brief outline of the methods will also be 
given. 

A correlated wave function for the neutron matter su- 
perfluid phase is constructed as 

I*,) = F |BCS) , (4) 

where the model BCS-state vector is 

|BCS> = Y[(u k +v k a^al kl )\0) . (5) 

k 

Uk and Vk are real, variational BCS amplitudes, satisfy- 
ing the relation + u| = 1, |0) is the vacuum state and 
a\ n is the fermion creation operator in the single-particle 
state, whose wave function is 

<p m =k,a(x = r,s) = —= r] a (s) exp(«fc • r) . (6) 

O is the normalization volume and 77 CT= |^(s) is the spin 
wave function with spin projection a. The second- 
quantized correlation operator F is written in terms of 
the N-particle correlation operators, Fjv, as 



F = E f» \K N )($ 



mjv I 
N I ' 



(7) 



where tun specifies a set of single-particle states. In 
coordinate representation and for a /6~type correlation, 
we have: 

( Xl ,X2,..X N \F N \$™ N ) = 

F 6 (l,2,..N){(f> mi (xx)<j> m2 (x2)..<f> mN (x N )} A . (8) 

The suffix A stands for an antisymmetrized product of 
single-particle wave functions. 

The cluster expansion of the two-body distribution 
function, (7(7*12), and of the one-body density matrix, 
n(rn'), for a simple Jastrow correlations (/,/(«, j) = 
f^( r *j)) was developed in Ref. 0. It is possible to 
sum, in a Fermi Hypernetted Chain (FHNC) theory, 
all the cluster diagrams contributing to the two quan- 
tities and constructed by the dynamical correlation lines 
(hj = fj — 1) and the by BCS statistical correlations, 

Po J (271T 
v f d?q 

= ft)7 (2^ ex P(* fe ' r ) u(r)u(r) , (10) 
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where po is the average density of the uncorrelated BCS 
model and v = 2 is the neutron matter spin degeneracy. 
For more complicated correlations, like the /g model, 
complete FHNC expansion cannot be derived. So, we 
have computed the expectation values at low orders of 
their cluster expansions. In particular, the average den- 
sity, 

_ (N) _ Em < a m a m > (u ) 

9 n n [ ' 

is computed at the first order of the expansion in a series 
of powers of the dynamical correlations. This expansion 
provides at each order the correct density normalization 
in the normal phase, since it fully takes into account can- 
cellations between same order diagrams but with differ- 
ent numbers and types of statistical correlations. Simi- 
lar, if not complete cancellations, hold in the superfluid 
phase. Consistently, the matrix elements of the Hamil- 
tonian on the correlated BCS state are evaluated at the 
two-body cluster level with vertex corrections at the in- 
teracting pair. 

In AFDMC for a BCS-like phase the neutrons are 
paired in a pfaffian 0> EH constructed by a trial BCS 
amplitude. To this aim, we have modified the nor- 
mal phase AFDMC code, substituting the Slater de- 
terminant with the pfaffian. We have taken the BCS 
amplitudes from the CBF calculations. The AFDMC 
simulation has been performed for N=12-18 neutrons 
in a periodic box and interacting through the Argonne 
vgi [21| (A8') potential. Preliminary results with larger 
N-values confirm the findings of this paper. Table (P) 
shows the results obtained by using a time step of At = 
5 x 10 -5 MeV -1 . The energy per particle for the normal 
phase at k-p = 0.6 fm for the closed shell case, N=14, is 
E NP /N = 2.548(3) MeV. For even N-values all neutrons 
are paired in the pfaffian, whereas for odd N-values the 
configuration of the unpaired neutron providing the best 
energy must be found. The gap energy for odd N-values 
is calculated according to: 

A(N) = E(N) - - (E(N + 1) + E(N - 1)) . (12) 

The present calculations show that, at shell closure 
for N=14, the normal phase is energetically slightly pre- 
ferred to the superfluid one. Two effects must be con- 
sidered: improvements to the pairing functions (e.g. a 
better choice of the nodal surface) and that the 14 par- 
ticles box dimensions are too small for the correlation 
length of the pairing function. Again, we have indica- 
tions that the superfluid phase is actually preferred from 
preliminary calculations with larger number of particles. 

The results for the gap function at the Fermi momen- 
tum are collected in Figure QJ. The curves labeled A° 
refer to the pure BCS estimates using operatorial reduc- 
tions of the Argonne vi% potential, down to the one. 



TABLE I: AFDMC energy per particle, E/N, and gap en- 
ergy, A, in MeV, in neutron matter with the ^48' potential 
for 12—18 BCS-paired neutrons at &f = 0.6 £m . In paren- 
theses the statistical errors are given. 



N 


E/N 


E 


A 


12 


2.6356(17) 


31.627(21) 




13 


2.7593(17) 


35.871(22) 


2.182(37) 


14 


2.5536(15) 


35.750(21) 




15 


2.8036(17) 


42.054(26) 


2.855(44) 


16 


2.6654(18) 


42.647(29) 




17 


2.8075(15) 


47.727(25) 


2.333(49) 


18 


2.6746(17) 


48.142(31) 





2.457(76) 



Most of the reduced potentials have been refitted to re- 
produce the S— and P-wave experimental phase shifts 
(V41-VS') ■ It is remarkable that all these potentials give 
the same BCS gap, provided they fit the same data set. 
Presumably this is a consequence of getting the effec- 
tive interaction near the Fermi surface correct. A ver- 
sion of Argonne vis simply cut to v§, without refitting, 
provides the higher BCS gaps of the A° 6 curve. In BCS, 
A° e = A° ig . The A„ 4 ,_ g , curves show the CBF gaps for 
the corresponding refitted potentials. The points with 
error bars are the AFDMC gap estimates with Argonne 
t>8'. The AFDMC results are the average values of the 
gaps 1|12[) computed around N=13, 15 and 17. 

The Figure shows that the low density (up to kp ~ 
0.6 fm" 1 ) CBF and AFDMC gaps, for the same A8' po- 
tential model, are in good agreement. The low order 
cluster expansion provides a CBF gap slightly smaller 
than the quantum MC one. The highest Fermi momen- 
tum AFDMC result, at kp = 0.8 fm" 1 , is much higher of 
the corresponding CBF gap, stressing the need of push- 
ing the cluster expansion to higher orders at increasing 
densities. It is clear, by looking at the A„ 4 , and A„ , , 
curves, that most of the reduction of the gap with re- 
spect to BCS is due to the tensor force, only partially 
compensated by the spin-orbit potential. The CBF A,j 6 
is the lowest curve, pointing once more to the importance 
of tensor force, in the sense that experimental data must 
be fitted in order to obtain meaningful results. 

The reduction of the gap due to medium effects is much 
smaller than previous estimates , which provided sup- 
pressions of a factor of two and more. The CBF theory 
shows an early disappearance of the gap, around p ~ 0.15 
Pnm, Pjva/=-16 fm" 3 being the empirical nuclear mat- 
ter saturation density. However, the QMC calculations 
do not seem to support this finding, and, surprisingly 
enough, give a gap energy that is still close to the stan- 
dard BCS one. A similar result was found in a prelim- 
inary study employing GFMC and simple model poten- 
tials H3|. 

The square in Figure gives the AFDMC gap, at kp = 
0.6 fm" 1 , for the A8' model implemented by the Ur- 
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bana IX three-nucleon interaction |2lJ. The gap with 
the three-body force results to be A Vgl+ mx—2. 810(146) 
MeV, slightly increased with respect to A„ , =2.457(76) 
MeV. A qualitatively different result was found in 
Ref. where the authors have obtained a small de- 
crease of the gap in a Bruckner G-matrix based ap- 
proach. 

In this letter we have used Quantum Monte Carlo and 
Correlated Basis Functions theories to microscopically 
evaluate the superfluid gap in pure neutron matter 
with modern interactions. These methods allow for tak- 
ing into account medium modification effects in a consis- 
tent and realistic way. In particular, QMC is expected to 
give a solution of the many-body Schrodinger equation 
very close to the true one Both theories are in good agree- 
ment up to the maximum gap density, and show a slight 
reduction of the maximum gap with respect to standard 
BCS. At higher densities QMC gives a larger gap than 
CBF, probably because the CBF gap is computed at low 
order of the cluster expansion and the missing diagrams 
become more and more relevant with the density. The 
gap reduction is essentially due to the tensor interac- 
tion. A novel and important effect is the small influence 
of the medium polarization, contrary to all the previ- 
ous estimates. The three-nucleon interaction provides a 
small increase of the gap. It will be most interesting to 
extend these analysis to other types of pairing, as the 
3 P2 — 3 F2 neutron pairing and the So proton pairing in 
highly asymmetrical nuclear matter. 



P [fm' 3 ] 

.00027 .00216 .00730 .01729 .03377 .05836 .09268 




FIG. 1: gaps for different nucleon-nucleon potentials and 
methods. See text. 
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